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ABSTRACT 

We use a semi-analytic circumstellar disk model that considers movement of the snow line through 
evolution of accretion and the central star to investigate how gas giant frequency changes with stellar 
mass. The snow line distance changes weakly with stellar mass; thus giant planets form over a wide 
range of spectral types. The probability that a given star has at least one gas giant increases linearly 
with stellar mass from O.4M to 3M . Stars more massive than 3M evolve quickly to the main- 
sequence, which pushes the snow line to 10-15 AU before protoplanets form and limits the range of 
disk masses that form giant planet cores. If the frequency of gas giants around solar-mass stars is 
6%, we predict occurrence rates of 1% for 0.4 M Q stars and 10% for 1.5 M© stars. This result is 
largely insensitive to our assumed model parameters. Finally, the movement of t he snow line as stars 
>2.5M© move to the main-sequence may allow the ocean planets suggested bv Le ger et al.l to form 
without migration. 

Subject headings: planetary systems: formation — planetary systems: protoplanetary disks — stars: 
evolution — stars: formation 



1. INTRODUCTION 

In the last ten years, the discovery of more than 200 
extra-solar planets, 2 and more than 200 debris disks, 3 
suggests that planet formation is a common and ro- 
bust process. Planet masses inferred from debris disks 
range from terrestrial to Jovian, at distances as great as 
tens of AU from the central star (e.g. lKenvon fe Bromlevl 
l2004al [Greaves et al.l 120051 ). The nature and sensitivity 
of radial velocity surveys means that most of the plan- 
ets are ~Jupiter mass gas giants in close orbits around 
Sun-like stars. However, recent discoveries as diverse 
as icy ^Neptune- mass planets orbiting M dwarfs (e.g. 
iRivera et~a l. 2005), and debris disks around A- type stars 
(e.g. lRieke et al.ll2005|) show that planet formation occurs 
over a wide range of spectral types. 

Current theory suggests that planets form in simi- 
lar ways around all stars. Thus, the increasing di- 
versity of stellar hosts and planetary systems provides 
an opportunity to test these theories. For this reason, 
the types of planets most likely to form around stars 
of differing s pectral types has become a renewed area 
of stu dy fe.g. Ilda fe Linl 120051 : lBossll2006t iKornet et all 
l2006t IKennedv et al.ll2006l ). after the ide a was first ex- 
plored by iNakanol nearly 20 years ago (|Nakanol 119871 
Il988allbh . 

Theories of Solar System formation generally include 
the "snow line," where ices condense from the nebular 
gas. The snow line distance is usually fixed in a disk 
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with a time independent surface density an d tempera- 
ture profile around a main-sequence star (e.g. Ilda fc Lira 
2005). In a more realistic picture, the disk and stel- 
lar properties evolve considerably during the 1-10 Myr 
pre-main-seque nce (PMS) lifetime when planets prob- 
ably form (e.g. iLissaueri fl987t iPollack et al.lll996| ). As 
the disk temperature evolves with time, movement of 
the snow line may therefore influe nce the properties of 
theoretical planetar y systems (e.g. IKennedv et al"1l2006t 
iGaraud fc Li^l2007l) . 

Here, we begin to develop a time dependent model for 
the formation of gas giant cores that considers the PMS 
evolution of the star and surrounding accretion disk. We 
introduce a simple semi-analytic disk model, based on 
the "minimum mass solar nebula," that links movement 
of the snow line through evolution of disk accretion and 
stellar luminosity. In contrast to pr evious studies (e.g. 
Ilda fc Linl 120051 : IKornet et al.l [2006) , our analysis sug- 
gests that gas giant formation around stars more mas- 
sive than the sun is more likely than around less massive 
stars. 

We cover the background important to our story in fj2j 
consider the snow line in <J2 and outline our model in <21 
We present our results in and discuss and conclude 
in Pand gZl 

2. BACKGROUND 

Planetary systems form in circumstellar disks, which 
evolve on timescales comparable to the pre-main- 
sequence (PMS) stage of stellar evolution. Ob- 
servations indicate a wide range of disk masses 
M disk ~ 01— 0.1 M+ (where M+ i s the stel- 
lar mass, e.g. lOsterloh fc Beckwithl Il995t iNatta et al. 
;n 2003 lAndrews fc Williams! 120051: lEisner fc Carpenter 
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[20061 iScholz et al.l 1200(1 a nd radii - 100-1000 AU 
(jMcCaughrean fc O The lifetime of the pri- 

mordial, optically thick, dusty component of the disk 
is 5,10 Myr, with a median timescale of ^3 Myr (e.g. 
IStrom et al.lll993t lHaisch et al.ll2001bl ). Though harder 
to observe, the gaseous compone nt of the disk is proba- 
bly re moved by viscous accretion I Lvnden-Bell fc Pringld 
19741) and photoeva p oration (e.g. Hollenbach et al.ll2000l : 
Adams et all [200l lAl cxander et alJ 120061)" on similar 



timescales ([Zuckerman et al. Hl995HPascucci et alj |2006) 
These timescales place strict observational limits on 
the important stages of planet formation. Planetesi- 
mals must form rapidly to enable further grai n growth 
and p rotoplanet formation by coagulation (e.g. ISafronovl 
1969). To attract significant atmospheres and form gas 
giant s, protoplanets need to reach masses of 5-10 M ffi 
(e.g. IPollack et all Il996t llkoma et all I2000D before the 
nebular gas is removed. 

In coagulation models, dust particles on near 
circular orbits with small relative velocities grow 
through repeated collisions and mergers in circumstel- 
lar dis ks. Further dynamical evolution through "run - 
away" (iWetherill fc Steward fl989l: iKokubo fc Idalfl99^ 
and "oligarchic" ( Kokubo fc Idalll998f ) growth leads to 
"isolated" protoplanets, whose mass Mi SO and spacing 
depend on their radial distance a from the central star 

via the Hill radius Rh = a (M isn / '!!,M i ,) 1 ^ (e.g. iLissauerl 
1987; Lis sauer fc Stevensodl2007f) 



M„ 



.3/2 



(3M*) 1 / 2 



(1) 



where a is the disk surface dens ity. Protoplanets ar e 
spaced at 2BR H ~ 8R H intervals (|Kokubo fc Idalll998D . 
Used in combination with equa tion (HI), the "niini- 
mum mass solar nebula" (MMSN, [WeidcnschilliH^ 11977] 



lHavashilH98ll ) with a 



oc a 



(where S = 1—1.5), gives a 



simple model of protoplanet formation. 

The "snow line" — the point in the disk that separates 
the inner region of rocky planet formation from the outer 
region of icy pla net formation — is an important feature of 
the MMSN (e.g ISasselov fc Lecarj [200q llda fc 03120051: 
iCiesla fc Cuzzill2006l) . Condensation of ices outside the 
snow line increases the disk surface density by a fac- 
tor fi ce ~ 3, 4 which leads to factor of 5 larger isolation 
masses (eq. [TJ. In an MMSN model with a = 10 gcm~ 2 
at 1 AU and 6 = 3/2, M lso w 0.1 (1) M e at 1 (5) AU. To 
achieve the probable cor e mass of 5-10 Mq for Jupiter 
(jSaumon fc Guillotl[2004T ) the MMSN can be augmented 
beyond the snow line by a factor ~4. Alternatively, if 
6 = 1 then M iso w 5 M m at 5 AU. Models that relax the 
assumption of a smooth radial profile fi nd an enhanced 
surfa c e density near the sn ow line (e.g. ICuzzi fc Zahnld 
2004: ICiesla fc Cuzzil 120061 A common theme among 
both MMSN and more detailed models is the surface 
density added by ice condensation. 

Because the timescale for planet growth is t oc P j a oc 
a 3 for a oc a " 3 / 2 , wh ere P is the orbital p eriod (e.g. 
ILissauerl [l987t see also iGoldreich et al.l (|2004f )). ice con- 

4 Th e usual value is ~4, but recent solar abundance figures for 
oxygen (Asplund ct al. 2005) indicate 3 is more reasonable. Recent 
composition d ata from 9P/Tempel f may argue for an even lower 
ice/rock ratio flKiippers et al.ll2005l ). 



densation also lea ds to shorter growth times . Numeri- 
cal simulations bv lKenvon fc Bromlevl (|2004al |bl) find the 
time to form 1000-3000 km objects agrees with this re- 
lation. Numerical estimates of the time to form the Jo- 
vian core range fr o m ~10 5 -10 6 yr (e .g. Lissauer 119871 : 
IPollack et all 11993 : llnaba et all 120031: IChambersll2T)03) . 
In general, the time to reach isolation ti SO provides an 
estimate of whether protoplanets form early enough to 
accrete gas and become giant planets. Short gas disk life - 
times (e.g. IZuckerman et al.lll995t iPascucci et alJfeOOSl) . 
imply a relatively short isolation time, and place strong 
constraints on the time to form gas giants by core accre- 
tion. 

Gas giant formation by core accretion occurs when pro- 
toplanet core masses are sufficient to attract gas from 
the nebula. T he core mass sets the timescale for gas g i- 
ant formation (Ikoma ct al. 2000; Hubick vi et al.l l2005). 
Cores with masses smaller than ~ 5M ffi attract atmo- 
spheres (e.g. llnaba fc Ikomal l2003f ). but are unable to 
form a gas giant before the nebular gas is removed on 
timescales of 1-10 Myr. Beyond the snow line the crit- 
ical core mass where significant gas accretion occurs is 
M crit ~ 7 M m M° f e K 25 (where M core is the rate at 
which planetesimals are accreted onto the core in units 
of 10~ 7 M m yr" 1 , and k is the gra in opacity in u nits of 
cm 2 g- MIkoma et al.ll2000L see also iRafikovl (|2006l )). The 
critical core mass re quired to form Jupiter in several Myr 



5-10 M a 



120051 ) — which implies n 



Pollack et ail Il996t iHubickvi etaLI 

10 g cm 2 at 5AU- 



-is consis- 
tent with th e core mass inferred from current structural 
models (e.g. ISaumon fc Guillotl[200l . 

2.1. Previous Work 

Most planet formation theories are based on a static 
MMSN disk around a solar-mass star. There are sev- 
eral motivating factors for extending these theories to a 
range of stellar masses: (i) the increasing stellar mass 
range of extra-solar planet hosts, (ii) observed trends 
with stellar mass, such as accretion rate and disk mass, 
and (iii) theoretical relations with variables that change 
with stellar mass, such as orbital period and isolation 
mass. This extension of Solar System theory to a range 
of spectral types began wi th a series of papers b y lNakanol 
nearly t wenty years ag o (iNakanol 119871 . Il9 8 8 al ibi). More 
recently, iKornet et al.l (|2006h considered formation of 
planets around stars of various masses in situ, while 
llda fc Lin (2005) examined observable planetary systems 
r esulting from t ype II migration. 

IKornet et alJ consider disk evolution prior to the 
growth of large objects. In their models, the increased 
inward migration rate for planetesimals around low-mass 
stars results in higher absolute surface densities from 0.1- 
100 AU at 1 Myr. Thus low-mass stars are more likely to 
form giant planets. This result is influenced by their 
choice of an approximately constant initial disk mass for 
all stellar masses. They do not consider planet formation 
b eyond 5 AU . 

llda fc Linl base their Monte-Carlo study on the MMSN. 
Type II migration — where a planet with sufficient mass 
opens a gap in the disk and whose orbit is subse- 
quently coupled to the v iscous evolution of the disk (e.g. 
iLin fc Papaloizoulll985l ) — is central to their model. In 
their attempt to reproduce the observed distribution of 
extra-solar planets, they find that close-in icy Neptune- 
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mass planets should be much more common than close- 
in Jupiter-mass plan ets around M dwarf s tars (see also 
lLaughlin et al.l [2004T )- In contrast to the iKornet et al.l 
study, they find that the likelihood of a star harbour- 
ing gas giants increases with stellar mass up to solar- 
mass stars. Their results are influenced by scaling the 
snow line distance as a snow cx M^, based on the main- 
sequence luminosity L+ cx M+. As we show below, this 
simplification places the snow line too close to (far from) 
the central star for stars with masses less than (greater 
than) a Solar mass when protoplanets form. 

In this paper we consider movement of the snow line as 
disk accretion subsides and the central star evolves to the 
main-sequence. Using our prescription for the snow line 
position over a range of stellar and disk masses, we locate 
regions where gas giant cores form. Assuming stars are 
born with disks from a distribution of masses, we then 
predict how gas giant frequency varies with stellar mass. 

3. LOCATION OF THE SNOW LINE 

In this section we consider evolution of the disk mid- 
plane temperature, and the snow line distance, with 
a simple model that includes accretion and PMS evo- 
lution. In particular, we are interested in the stellar 
mass dependence, rather than a detailed derivation for 
a single star. As we show in <j5j ^1AU differences 
between our model and more detailed tr eatments (e.g. 
ISasselov fe Lecarl 120001: iLecar et al.ll2006D do not affect 
our conclusions. 

The disk mid-plane, where the gas density is highest, 
is probably where most ices condense, and has a temper- 
ature 

+ Tf rr , (2) 



rp4 

mid 



rrii 

mid,accr 



where T m id. a ccr is the mid-plane temperature arising 
from viscous forces within the disk, and Ti rr is tempera- 
ture due to external irradiation of the disk by the central 
star. 

The effective disk temperatur e from viscous accretion 
is (|Lvnden-Bell fc Pringldll97l 



T. 



3 GM+M 



eff,accr 




(3) 



where M is the accretion rate, is the stellar ra- 
dius, and <j s b is Stefan's constant. In optically thick 
regions near the snow line, the mid-plane t empera- 



turc is T„ 



mid.accr 



3 T T e ff,accr/ 



(|Hubeny|[i990l ). where 
r = na g 1 2. The opa city k is a function of tempera- 
ture (jBell fc Linlll994f ). and the gas surface density a g 
is 100 times greater than that of solids. The figure of 
100 is used in conver ting mm dust obse rvations to to- 
tal disk masses (e.g. iNatta et al.l I2000D based on the 
interstellar gas/dust ratio, and i s similar to the sola r 
metallicity (Z) fraction of 0.0122 (|Asplund et al.ll2005[ ). 
The accretion rate varies with stellar mass as approxi- 
mately M cx M+ for the range of stella r masses we con- 
sider (0.2-4 M*, iMuzerolle etaH l2003f) . and with time 

as M cx (t/10 6 yr)" 7 . lHartmann et alj (| 19981 ) derive 
7 = 1.5-2.8. The uncertainty is due to the limited age 
range of their sample and a large range of accretion rates 
at a given age. The value 7 = 1.5 is their "preferred re- 
sult." We scale M with surface density, which accounts 



for the observed trend with stellar mass (if disk mass 
scales linearly with stellar mass, see @, and is consistent 
with expected viscous evolution (where M cx va and v 
is the disk viscosity). For ~lMyr old Solar- type stars 
M ~ 10- 8 M Q yr" 1 (jHartmann et al.l fl99l . We set 
M = 10 -8 Mq yr -1 for an initially three-fold enhanced 
MMSN disk, as this disk decays to the "typical" observe d 
~MMSN mass disk by ~lMyr (jHartmann et al.lll998[) . 

A more complete treatment of the optical depth to the 
mid-plane would include evolution of the gas surface den- 
sity, allowing the mid-plane temperature to drop some- 
what faster than described above as a a decreases in the 
inner disk (Lvndcn-Bel l fc Pringldll974l ). The solid sur- 
face density, which largely resides near the mid-plane and 
determines protoplanet characteristics, remains largely 
unaffected by the gas disk evolution (aside from snow 
line evolution). 

The disk temperature contribution from irradiation is 



a\ 1/4 (R, 



3/4 



(4) 



where a « 0.005/dAU + 0.05a A jJ fo r a flared disk in ver- 
tical hydrostatic equ ilibriu m (e.g. lAdams fc Shu l fl986h 
iKenvon fe Hartmannl [l987l IChiang fc Goldreichl 119970 . 
Here oau is a in units of AU. When the disk is opti- 
cally thick to radiation at this temperature, T irr is ap- 
proximately the interior t emperature for a flared disk 
(|Chiang fc GoldreichlH997l) . 

Figure [T] shows the location of the snow line in disks 
with a = ctmmsn M*/Mq for several different stellar 
masses over time for irradiation only, and for accretion 
+ irradiation. We locate the snow line where T m id = 
170 K. More detailed derivations of this tem perature (e.g. 
iPodolak fc Zuckerll2004l ; ILecar et al]|2006f ) do not change 
the snow line dist ance significantl y . Fo r PMS stellar 
properties we use iPalla fc Stahlerl (|1999( ) tracks. For 
comparison, we also show the (fixed ) snow line dista nce 
for stars on the main-sequence (e.g. llda fc Lhl[2005l) . 

Looking first at the Solar case, the snow line moves in- 
ward over time. This movement is always determined by 
viscous accretion, and its decay over time. Our snow line 
crosses the "canonical" distance of 2.7 AU at 5 x 10 5 yr. 
For disks with accretion rates so low that irradiation 
dominates, the snow line still moves inward over time, 
as illustrated by the dashed line for Ti rr . 

For more massive stars, Tj rr begins to dominate as 
accretion subsides and the star quickly evolves to a sig- 
nificantly greater main-sequence luminosity. Irradiation 
becomes important at a few Myr for 2 M stars, and 
~lMyr for 3M stars. T he large discrepan cy between 
our snow line, and that of llda fc Lin! l|2005l ) (who con- 
sidered 0.2 M Q < M+ < 1.5 Mq) arises because theirs is 
based on the main-sequence luminosity (L+ cx M%) and 
an optically thin disk {Tj isk cx L* a~ 2 ). With our model, 
the snow line distance is less sensitive to stellar mass, 
allowing icy protoplanet formation relatively close (~5- 
10 AU) to the central star for intermediate mass stars. 
At these closer distances, the surface density is higher 
and formation is faster, making it more likely that pro- 
toplanets massive enough to undergo core accretion will 
form. For less massive stars, the snow line is still at a few 
AU, where isolation times are relatively long, making it 
difficult to form cores before the gas disk is dissipated. 
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Fig. 1. — Location of the snow li ne (a at T m id = 170 K) over time for 0.6, 1, 2, and 3 Mq stars (left to right, and down) with irradiation 
only (using Palla & Stahler (1999) PMS tracks, dashed line), and irradiation + accretion (solid l ine). The disks have surface densities 
o" = ctmmsn^*/^©- Included for reference IS CLsnow — 

2.7M*/M AU as used bv llda fc Linl (2005) (dotted line). 



Comparison of the snow line distance with typical disk 
lifetimes of several Myr leads to an increasing snow line 
distance with stellar mass. 

With the snow line evolution established, we now de- 
scribe our model of protoplanet formation. 

4. PROTOPLANET FORMATION MODEL 

The MMSN is a simple model disk for the origin of 
the Solar System, and has a (a) = <J n f ice a~ s , where 
the factor fi ce represents a jump in surface density at 
the snow line distance a snow , and S is usually 3/2. To 
extend this model to a range of stellar masses requires 
consideration of how disk mass varies with stellar mass . 
Observations indi cate Mdisk oc Ah (Natta et al.1 [2000; 
iScholz et~a l. 2006); however, there is a wide range of disk 
masses at any given stellar mass. Thus, to extend the 
MMSN model to a range of stellar and disk masses, we 
adopt the surface density relation 

M* , 



a{a, t) = a t) fu 



(5) 



where <tq — 10 g cm~ 2 . The factor rj changes the disk 
mass relative to the star ("relative disk mass"), and is 
varied to account for the observed range of disk masses at 
fixed stellar mass. Current observations suggest rj ~ 0.5- 
5 (M disk = 0.01-0.1 M*); r] ~ 10 is the upper limit for 
disk stability {Mdisk ~ 0.25 M*). To provide a smooth 
transition from /j ce = 1 for a < a snow to /, ce = 3 for 
a > a snow , we set f ice = 1 + (A,; ce - 1)/(1 + e x ) where 
^■■ice — 3, x — {oj Snow cl) I ' Aci snow and Aa sf)ot[J is the 
radial distance equivalent to a 5 K temperature change. 



Combined with the local orbital period, the surface 
density sets the time to form protoplanets, and sets 
whether protoplanets form early enough to accrete gas 
and become gas giants. We introduce a stellar mass de- 
pendence, so o ur isolation timescal e . based on numerical 
simulations bv lKenvon fc Bromlevi (|2004al rbl) . becomes 

t iso oc (wy 1 a 3 / 2 A/7 1/2 . (6) 

The normalisation of e quation §E§ depends on the size of 
the s mall objects (e.g. iGoldreich et al.ll2004t IChambersl 
120061 ). We use 10 5 yr for a = 10 g cm" 2 at 5 AU, 
based on the likelihood of small fragm ented bodies to 
accrete (e.g. iKenvon fc Bromlevi l2004al). and the con- 
seque nt short growth t imes ( Rafikovl 12004 IChambersl 
l200l . iHubickvi et all (f20Q5h infer t iso <5xl0 5 yr for 
much larger 100 km planetesimals. As long as the Jovian 
timescale is somewhat shorter than the gas disk lifetime, 
this choice affects our results little. 

Under the assumption that stars all form disks 
in ^1 5 yr (based on an infall rate of 1O -5 M0 
yr -1 , iPalla fc Stahler! I1999T ) . and that plane t esimal 

formation is relative ly fast ( Wcidcnschil lind 120001 : 

iDullemond fc Dominikl (20051) . we add a constant offset 
of 10 5 yr to the isolation timescale to reconcile the tim- 
ing of isolation with disk and stellar evolution. Though 
this time is uncertain, removing or moderately modifying 
the offset does not affect our results significantly because 
ti SO is usually >10 yr. 

We adopt a range of disk masses, integrated from the 
inner disk radius to 60 AU, with a gas to solids ratio 
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of 100. For S = 3/2, rj — 4 corresponds to a rel- 
atively massive disk Mdisk — 0.1 Af*. This enhance- 
ment is our baseline model and yields the surface den- 
sity and cor e mass needed to form Jupiter on reaso nable 
time scales ([Pollack et al.lll996t Hkoma et alJl200Ct ). For 
5=1, smaller rj yields the same disk mass because more 
mass is placed at larger radii. For rj = 1 and 5 = 1, 
M^sk = 0.12 Af*. 

5. REGIONS THAT FORM GAS GIANT CORES 

The successful formation of a gas giant planet by core 
accretion requires satisfaction of two main conditions. 
A core must form while the gaseous component of the 
circumstellar disk is still present, and it must be massive 
enough to attract a large atmosphere before this gas is 
dispersed. Prior to isolation, accreted planetesimals and 
a sub-critical protoplanet mass limit gas accretion. After 
isolation, if a protoplanet is massive enough, and forms 
while the gas disk is still present, significant gas accretion 
proceeds. This separation into two classes, gas giants 
and "failed cores," reflects the expected pau city of 20- 
100 M m planets over a range of stellar masses (llda fc LinI 
[2001 . 

To form a gas giant in less than 10 7 yr, va rious studies 
suggest a minimum core mass of 5-10 Afcp (|Ikoma et all 
l2000Hlnaba et al.ll2003t iHubickvi et al.ll2005f ). Although 
the core mass depends on the planetesimal accretion rate 
and the opacity, the derived sensitivity is weak (M core oc 
AT - 25 ; llkoma et al.ll2000L see also iRafikovl ([200l i A 
limited reservoir of planetesimals to accrete after isola- 
tion means high accretion rates cannot be sustained (and 
M will decrease), while low accretion rates lower M core . 
Thus, we adopt a minimum core mass M core = 10 Af® as 
a baseline, and consider M core = 5 Mq in £15.41 

The timescale for gas dissipation sets our second re- 
striction. The gaseous component of the disk disperses i n 
<10Myr (jZuckerman et all 119951 : iPascucci et al.ll2006f ). 
With Af oc a, the dissipation timescale for a viscous disk 
td oc Mdisk/ M ~ constant for our assumptions. Be- 
cause the disk mass decreases significantly (~60%) in 
1 Myr, we adopt t core = 1 Myr as a typical maximum 
core formation time for all disks. Henceforth we reserve 
the word "core" for a protoplanet with Mj so > M core and 
tiso < tcore.- Relatively little is known about the evolu- 
tion of the gaseous component of the disk; we comment 
further on the consequences of varying t core and other 
parameters in Q5A\ 

To investigate locations within circumstellar disks 
where gas giant cores form, we first derive results for 
a 1 M Q star, and then consider a range of stellar masses. 
We restrict our study to stars with masses 0.2-4 Mq. 
For stars with masses <0.2 Af© our model does not form 
gas giants. The short main-sequence lifetime of massive 
stars (Af* > 4M a ) makes them much less likely tar- 
gets for planet detection. Some oligarchs do not reach 
masses sufficient for core accretion. Those cores that do 
form compete with other cores for dynamical space in 
the disk. We defer consideration of these objects to S 35.3I 

5.1. The Solar Example 

In this subsection we show how the MMSN disk model, 
the moving snow line, and the isolation mass and time 
combine to give a picture of the Solar System structure 
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Fig. 2. — Isolation mass (filled circles, labelled with Mi ao ) as a 
function of radial distance and PMS model time, for a solar mass 
star with the MMSN model with rj = 4 and 8 = 3/2. Masses are 
spaced at 8 Rh intervals and only shown outside the snow line. 
The solid line shows a s ?iow over time, and the dashed vertical line 
is the time t core = 1 Myr. 

at ~lMyr. 

Figure [2] shows isolation masses for the MMSN model 
beyond the snow line with rj = 4, as a function of time 
and radial distance from the Sun. The isolation mass and 
timescale are calculated from the equ ations described 
in $S Accretion and PMS tracks from iPalla fc Startler! 
(1999) set a snow as described in Sj3l 

When the first objects reach isolation, the Sun is in 
the early stages of its PMS contraction, and the accre- 
tion rate is ~10~ 7 Mq yr -1 . Consequently the snow line 
is at ~6 AU, which determines where the innermost icy 
protoplanet forms. This protoplanet is massive enough 
to become a gas giant (a core) , so we refer to this position 
as the inner edge of the core-forming region. In the ab- 
sence of significant migration from disk interaction, this 
result may help explain why Jupiter is at 5 AU. As the 
Sun continues to contract and accretion decreases, isola- 
tion is reached at ever increasing distances beyond the 
snow line. Eventually, the isolation time becomes longer 
than t core , and protoplanets form too late to undergo 
core accretion. Isolation masses increase with distance 
from the Sun, so t core always sets the outer edge of the 
core-forming region, and the number of cores that form. 
If the cores are spaced by 8 Rh at isolation (as in Fig- 
ure [2]) , then ~4 cores form in this region and the region 
extends from ~6-ll AU, similar to the region containing 
Jupiter and Saturn today. 

5.2. A Range of Stellar Masses 

We now consider how the core-forming region for a 
solar mass star changes with stellar mass between 0.2 Mq 
and 4Af©. The processes described in S}5T] still apply, 
but differences arise due to the linear relation between 
disk surface density and stellar mass, the different orbital 
periods around other stars, and the changing snow line 
distance due to the evolution of accretion and the central 
star. 

Figure [3] shows isolation masses for 0.6 and 3 Mq stars, 
with t core = 1 Myr and the same relative disk mass as 
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Fig. 3. — Same as Figure [2] but for 0.6 Mq (left) and 3Mq (right). Isolation masses are only plotted outside the snow line, or where 
Mi SO > M core . 



Figure [2j The lower mass star does not form any cores, 
as the first object with M iso > M core forms after t core . 
However, large 5-10 Mq objects still form (see ^5.3.1|) . 
A longer t core allows some of these to become cores, so 
whether lower mass stars form gas giants is sensitive to 
both Mcore and t core . 

The 3 Mq star forms its innermost core just inside the 
snow line at ^8AU, and the outermost is at ~20AU. 
Cores can form before t core at greater distances due to 
decreased P and increased a. The greater surface den- 
sity in disks around these stars all ows rocky cores t o form 
interior to the snow line (see also llda fc Linll2005f ) . The 
large mass of the cores (^>M core ) in these relatively mas- 
sive disks probably allows some gas accretion prior to 
isolation, which requires a numerical model of growth 
for more investigation (see also §5.4[) . 

The results over a range of stellar masses can be com- 
bined into a single figure that considers the core-forming 
regions as a function of stellar mass. Figure |4] shows the 
core-forming region as a function of radial distance and 
stellar mass for our model, with standard parameters of 
tcore = lMyr, Mcore = 10 M®, and S = 3/2. Each con- 
tour represents the inner, outer, and stellar mass limits 
for forming cores with a particular relative disk mass. On 
the r\ — 4 contour (outlined in the figure), the 6-10 AU 
range from Figure [2] contributes the points y = 6 AU and 
y = 10 AU for x = 1 Mq. Similarly, the 8-20 AU range 
from the right panel of Figure [3] contributes points at 
x = 3M Q . 

In general, Figure |4] shows that as stellar mass in- 
creases, disks with lower relative disk masses form cores, 
and the width of the regions where these cores form 
increases. Stars more massive than ~12M form 
rocky cores interior to the snow line for relatively high 
disk masses. The core-forming region expands outward 
with increasing relative disk mass because the isolation 
timescale becomes shorter. Doubling the disk mass al- 
lows a to increase by a factor of ~1.6 to keep the same 
Uso (eq. EJ. 

As stellar mass increases, cores form in disks with de- 
creasing relative disk mass. As rj decreases, the inner 



edge of the core-forming region moves inward, because 
the accretion rate is lower and the disk has lower opti- 
cal depth (so the mid-plane is cooler), and the snow line 
has evolved closer to the star by the time isolation is 
reached. The snow line is roughly the lower edge of the 
darker (lower rj) contours. 

If Mi SO does not jump to a value >M core due to /,* ce 
at the snow line, a core still forms further out. Thus, 
the inner edge moves to greater distances as the stellar 
mass — and hence absolute disk mass — decreases for fixed 
relative disk mass. However, the t CO re restriction means 
that oligarchs around sufficiently low-mass stars do not 
reach isolation in time, and that all contours have a lower 
stellar mass limit. 

The lowest disk mass that forms cores has rj = 1.2, and 
only does so for 3M . For more massive stars, irradi- 
ation overcomes accretion as the star reaches the main- 
sequence, and the larger snow line di stance makes for - 
mation of cores more difficult (see also llda fc Lin1l2005h . 
Thus 3 Mq stars are the most likely to form at least one 
gas giant, as they form cores over the widest range of 
disk masses. 

The width of the regions w here cores f orm i ncreases 
with stellar mass (Figure H llda fc Lin! I2005D . The 
spacing of cores remains roughly constant with differ- 



ent stellar mass (Mj. 



OC (7 



AM* oc M* and Rh oc 



(Mi so /M*)V 3 for 5 = 3/2, see also Figs. [2]&[3|. How- 
ever, the number of cores is not linearly related to the 
region width, since the spacing becomes wider with in- 
creasing distance. The width of the regions depends 
strongly on disk mass, particularly for disks that form 
cores interior to the snow line. The increasing width 
of the core-forming regions suggests that the number of 
cores (and therefore planets) in individual planetary sys- 
tems increases with stellar mass. 

To summarise, the range of relative disk masses that 
form gas giant cores increases with stellar mass, as does 
the width of the regions they form in. The first result 
leads to the expectation that the likelihood of forming 
gas giants increases with stellar mass. We make a quan- 
titative prediction in £16.11 
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While the core-forming regions are our primary inter- 
est, there are large regions of parameter space where oli- 
garchs are relatively massive, but will not form gas gi- 
ants. We consider these planets now. 

5.3. Other Planets 

The discovery of exoplanets with masses smaller than 
Neptune suggests that planet formation might often yield 
failed cores — oligarchs that did not accrete gas from the 
disk. Neptune and Uranus may be co nsidered failed 
cores . Several theore tical studies (e.g. lLaughlin et al.1 
120041: llda fc Lin! 120051 this paper) suggest that failed 
cores are more common around low mass stars. 

5.3.1. Failed Cores 

In our model, there are three ways to produce failed 
cores. Objects that form too late (ti SO > t core ), or form 
with insufficient mass to accrete gas (Mi so < M core ) are 
failed cores. Although all cores within the core-forming 
region can potentially accrete gas, dynamical interactions 
among the cores may eject one or more into regions with 
a small gas surface density (e.g. iThommes et al.l Il99 
or from the system entirely (e.g. iGoldreich et al 1 1200 
iFord fc Chi"angll2007f l This mechanism occurs m a ran- 
dom (and currently unquantifiable) fraction of models 
with Mi SO > M core and ti SO < t core . In our model 
failed cores are more common around low-mass stars be- 
cause isolation masses are smaller, and isolation times 
are longer. 



Apparent fail ed cores may also for m by collisions over 
long timescales (jKennedv et al.l l2006). The mass of these 
icy planets may be limited by the likelihood of colli- 
sions vs. ejections d uring the final stages of coalescence 
(|Goldreich et al.ll200l . 

5.3.2. Ocean Planets 

The dive r sity of observed extra-solar planets led 
ILeger et"aLl (|2004l) to suggest that 1-10 M® icy planets 
that form in the region beyond the snow line may migrate 
inward to ~1AU, where the outer layers subsequently 
"melt." With masses too low to accrete much gas, these 
planets are less dense than a rocky pla net of equiva lent 
mass and harbour deep oceans. Hence ILeger et al.l call 
these "ocean planets." 

The increase in luminosity of stars with masses 
>2.5Mq as they reach the main-sequence provides an 
alternative in situ formation mechanism for ocean plan- 
ets, as they may have insufficient mass for significant 
migration. We outline the concept briefly, because these 
planets are difficult to detect. At times >l-10Myr, the 
snow line moves to ~10 AU as the disk becomes optically 
thin. 5 Failed cores in the range 1-10 M® can therefore 
achieve their final mass outside the snow line in ~1 Myr, 
and without migrating, later find themselves in a much 
warmer region when the star reaches the main-sequence. 

5 The snow line has less meaning at these times, since there is 
little gas to condense into ices. The equilibrium temperature of 
objects is a more relevant concept. 
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Fig. 5. — Same as Figure [4] but for 8 = 1 and M core = 5M . 
For r) = 0.25, only 3 Mq stars form cores. Though the <5 = 1 disk 
has different r), it covers the same range of disk masses as Figure [4] 
with 8 = 3/2, with the exception of the lowest disk mass (rj = 1.2) 
from that figure. 

Though all stars more massive than the Sun undergo 
an increase in luminosity as they set tle onto the main- 
sequence (e.g. iPalla fc Stahlerl Il999). the temperature 
at the early snow line distance of ~7AU must increase 
enough to melt ice and maintain oceans. For 2.5 Mq, 
there is ~1 AU overlap between the early snow line, and 
the final habitable z one distance — at an equilibrium tem- 
perature of ^245 K (Kastin g et al.lfl993h — with room for 
a few cores that form in situ just beyond the snow line. 

5.4. Sensitivity to Model Assumptions 

Our model is simplified, but captures some important 
concepts. In this section, we show that our results remain 
for realistic variations on our model assumptions. 

We use a simple model for the temperature profile of an 
irradiated, accreting, flared disk, which sets the location 
of the snow line. As shown in Figured! ~AU changes in 
the snow line distance affects where the innermost cores 
originate, but there is little change in the range of disk 
masses that forms massive cores for a given stellar mass. 

The disk surface density profile is uncertain: the 
MMSN assumes the Solar System planets formed in situ. 
In the standard MMSN model, (5 = 3/2 and rj = 4 yield 
the surface density needed to form a massive core near 
Jupiter. If 8 = 1 and n = 1 , then M iso 5 (13) Mq at 
5 (10) AU, which are si milar to the inferred cor e masses 
for Jupiter and Saturn (|Saumon fc Guillotll200l . Fi gurc 
[5] shows the core-forming regions for T) = 1 , 8 = 1 and 
M core = 5 Mq . The lower a in the core-forming regions 
needed to keep the same disk mass makes it harder to 
form cores with M core = 10 M®, but with M core = 5 Mq 
the regions are similar to our baseline model. 

There is little observational constraint of gaseous inner 
disk lifetimes, so the least certain of the parameters we 
specify is t core . The isolation time is ti SO oc a? (when 
8 = 3/2), so doubling t core allows the outer edge to move 
outward by a factor of about 1.3. Figure [6] shows how 
changing t core alters the core-forming region with r\ = 4 
for M core — 5 Mq and 10 M®, and 8—1 and 3/2. Longer 
gaseous disk lifetimes lead to more gas giant cores. The 
general trend is to extend the regions to lower stellar 
masses and to greater radial distances. Cores that take 
longer to form — due to smaller P or a — can reach isola- 
tion before the gas disk is dissipated. 



Dust disks around spectral types earlier than ^G have 
some what shorter lifetimes than disks around lower mass 
stars (|Haisch et al.ll2001al ICurrie et al1l2007t ). If gas is 
removed on timescales similar to infra-red excesses for 
a range of spectral types, t core is shorter and the outer 
edge of the core-forming region for these stars moves in 
(Fig. [S]). Observations of evolved stars can test whether 
a strongly stellar mass dependent t core (or some other 
process) defines an upper stellar mass limit for gas giant 
formation. 

The short timescale for disk removal by the cen- 
tral s tar — inferred from the lack of transition disks 
(e.g. iKenvon fc Hartmlnml 119951 : IClarke et all 1200 It 
lAlexander et alj 12006ft — and the likely large radial dis- 
tance in the disk fo r external influence by massive stars 
([Adams et al.l [20041. supports our assumption that the 
location of photoevaporation is unimportant. Although 
the expected distance for external photoevaporation in 
"typical" cluste rs reaches closer to the central star for 
low- mass stars (I Adams et al.ll2003 ). it lies outside the 
core-forming region for our standard case of M core = 
10 M®, t core — 10 6 yr, and 8 = 1.5. The core-forming 
and external photoevaporation regions begin to overlap 
for M* < 0.5 Mq and t core > 5 x 10 6 yr. 

The size of planetesimals is uncertain, as is the iso- 
lation time that results from their accretion by oli- 
garchs. Our choice of 10 5 yr for the Jovian core is 
relati vely short, and is based on likely fragmentation 
(e.g. IKenvon fc Bro mley 2004ah and the rapid accretion 
of small ~1 00m planetes i mals in the she ar dominated 
regime (e.g. lRafikovll2004t IChambersl llOOe). If planetes- 
imals are larger and t{ so is longer (e.g. lMyr), similar 
results can be obtained by simply using a longer (yet still 
reasonable) t core ~ 3 Myr. 

Choosing M core = 5 Mq allows core formation in less 
massive disks. Halving M core allows a disk 1.6 times 
(Mi SO cx <7 3 / 2 ) less massive to form cores in the same 
region, and also extends the region to lower stellar masses 
(Fig.®. 

We have chosen to ignore type I migration, where lin- 
ear theory predicts that protoplanets excite spiral density 
waves in the gas disk, and m igrate toward the ce ntral star 
on 0.01-0.1 Myr timescales (|Tanaka et al.ll2002f) . Recent 
studies indicate that for cores less massive th an ~10M® 
the timescale is longer (jMasset et al.l 12006) , and may 
be reduced to a r andom walk due to magne tohydrody- 
namic turbulence ([Nelson fc Pap aloizou 2004) . For cores 
with masses ^10 Mq (Fig. [3|) in relatively massive disks 
around intermediate mass stars, core accretion will likely 
occur bef ore isolation, w hile planetesimals arc still being 
accreted ([Rafikovl 12006). The successful formation of a 
gas giant then depends on whether the planet can reach a 
gap-opening mass before migrating into the central star. 

Once a planet opens a gap in the disk it has survived 
the type I migration regime, but its continued existence 
is not guaranteed. Depending on the disk viscosity and 
lifetime, the planet can still migrate onto the central star 
by type II migration. 

In summary, the simplicity of our model means that 
reasonable changes in the input parameters change the 
results little. Future development of the model can in- 
clude a more complete treatment of more complicated 
physical processes. 
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Fig. 6. — Similar to Figure [4] but for fixed r] = 4. Contours represent different t core as indicated by legends. Top (bottom) panels are 
6 = 1.5 (<5 = 1), and left (right) panels are Af core = 10 (M core = 5Mj). 

6. DISCUSSION 



The age of direct planet detection is approaching (e.g. 
NICI Campaign on Gemini South) , where discoveries will 
be pushed to larger semi-major axes. In addition, ra- 
dial velocity surve ys now extend over a wider range of 
stellar masses (e.g. iFrink et al.ll2002t iButler et al.ll2006t 
I Johnson et all l2007al jbl).~ Our goal is to develop a the- 
ory of planet formation that extends over the observa- 
tional range to make testable predictions, and to develop 
greater insight into the processes that produce the ob- 
served diversity of planetary systems. 

For planets orbiting giant stars, there is a downward 
shift in the planet-m etallicity distribution by ~0.3 dex 
(jPasquini et al.l l2007f ). This result is not surprising in 
the context of our model. In Figure |3J the lowest rel- 
ative disk mass that forms cores roughly halves from 
1 to 2Mq, which corresponds to a —0.3 dex change in 
metallicity. Thus, we naturally expect the lower end of 
the metallicity distribution of higher mass stars to be 
shifted. However, we do not expect the high metallicity 
end of the distribution to move, since these disks can still 
form cores. 

There are now sufficient planet discoveries to start 
quantifying trends across a range of stellar masses, which 
allows the first steps towards comparison with planet 
formation th eories that cons ider the mass of the cen- 
tral star (e.g. Ilda fc Linll2005| ). Though sample numbers 
are small, studies of >1.3M Q giants indicate that giant 
planet frequency increases with stellar mass in the range 



0.1-2 Mq ([Johnson et al.ll2007af ). We now calculate what 
our model predicts for the probability of forming gas gi- 
ants as a function of stellar mass. 

6.1. Gas Giant Frequency and Stellar Mass 

Assuming all stars are born with a distribution of disk 
masses, we can estimate the probability Pqg of a star 
forming at least one gas giant as a function of stellar 
mass. Though compari son with observe d disk masses 
is uncertain, we follow Ilda fc Lin] (|2005f ) and adopt a 
Gaussian distribution in terms of x = log Mjisk/M* , 

where Pdisk oc exp [x — fi) 2 /2<rf n ^ with standard de- 
viation cri n = 1/3, centered on Aldisk — 0.03 M* (e.g. 
[i k, —1 .5). This d i stribu tion is similar to data com- 
piled bv lNatta et al.l (|2000f ). which is sensitive to all disk 
masses that form a core in our baseline model. Because 
our model is based on parameters that change with stel- 
lar mass (such as isolation time and disk mass), the rel- 
ative probability of forming gas giants is our main con- 
cern. Effects that may set the absolute probability, such 
as survival of migrating planets, are not included. To 
make contact with observations, we t herefore normalis e 
our results to 6% for Solar-mass stars (|Udrv et al.l [2007) . 

Figure [7] shows the likelihood of a star harbouring at 
least one gas giant planet as a function of stellar mass for 
our baseline model (FigH|) and a model with (5 = 1 and 
M core = 5M ffi (FigEJ), normalised to 6% at 1 M . To il- 
lustrate the difference between a static, main-sequence 
scaled snow line, and our evolving one, we include a 
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Fig. 7. — Probability of a star harbouring at least one gas giant 
planet as a function of stellar mass for our baseline model (solid 
line), and <5 = 1 and M core = 5M(j (dashed line). The thin 
dotted line is a fitted line of constant slope Pm« = 0.20 M* — 0.06. 
The dot-d a shed line has a anow oc 2.7M^ AU, for comparison with 
llda fc Linl 12005D . All curves are normalised to 6% at 1 Mq via a 
straight line fit. 

model with a S now oc 2.1 Ml AU, similar to the main 
model of llda fc Linl (|2005| ). Each point is the probability 
a star has a disk in the range that forms cores in Fig- 
ured] This plot assumes that if one or more cores form, 
at least one will result in a gas giant. For comparison to 
current observations, it also assumes that as they accrete 
gas, cores generally migrate or scatter a stellar-mass in- 
dependent fraction to observable distances, where the 6% 
normalisation applies. Like our baseline model, different 
M core , t core , and 5 have lines of approximately constant 
slope {Pgg = mM* — c) up to ~3M Q . Our baseline 
model has m = 0.20 and c = 0.06; very few gas giants 
form by core accretion below 0.3 Mq. 

With the normalisation, our baseline model predicts 
1% of 0.4 Mq stars, and 10% of 1.5 Mq stars will har- 
bour at least one gas giant. Increasing ai n decreases the 
range of probabilities, because disk masses come from a 
less strongly varying section of the overall distribution. 
For example, cri„ — 1 yields 4% and 8% for 0.4 M Q and 
1.5 M stars respectively. 

Our result is robust to changes in our assumed model 
parameters (in the range 0.4-1. 5M©), because the range 
of core forming disk masses generally remains the same, 
and the 6% normalisation removes absolute differences 
for different model parameters. 

Despite its simplicity, our model modified to have a 
constant snow line at a snow = 2. 7 (M+/ M P) ) 2 AU pro- 
duces a similar result to llda fc Linl {2005) : that less gas 
giants form above a Solar mass. The difference arises 
from the stronger dependence of snow line distance with 
stellar mass. For stars <3M in our model, accretion 
largely determines the snow line distance, which sug- 
gests a better scaling is a snow oc M* a^ S if M oc <tq and 
o g oc A/*a~ 3 / 2 . Alternatively, a snow oc M^ 9 if M oc M+, 

or a snow oc M* /9 if M oc M*. 

Finally, our probability calculation does not take the 
increasing core-forming region width with increasing stel- 
lar mass into account. Cores that form later at larger 
distances may be less susceptible to migrating into the 



central star, and stars with multiple planets may be more 
likely to retain at least one during migration and scatter- 
ing processes. These effects have the potential to increase 
the frequency of giant planets as stellar mass increases. 

7. SUMMARY AND CONCLUSIONS 

We describe a model for the evolution of the snow line 
in a planet forming disk, and apply it over a range of 
stellar masses to derive the probability distribution of 
gas giants as a function of stellar mass. The two main 
ingredients for our model are a prescription for movement 
of the snow line due to accretion and PMS evolution, and 
rules that determine whether protoplanets are massive 
enough, and form early enough, to become gas giants. 

The snow line distance generally moves inward over 
time. With our prescription for the accretion rate, ac- 
cretion dominates over irradiation for stars with M* < 
2Mq. For >3Mq stars, irradiation dominates at times 
>1 Myr as the star moves up to its main-sequence lumi- 
nosity. The transition is at a few Myr for ~2 Mq stars. 
Over the wide range of observed accretion rates for any 
fixed stellar mass, the snow line in some disks may be set 
entirely by irradiation. 

The snow line generally sets where the innermost gas 
giant cores form. In relatively massive disks around in- 
termediate mass stars, rocky cores form interior to the 
snow line. The location of the outermost core is always 
set by the gas dissipation timescale. The range of disk 
masses that form cores, and the radial width of the re- 
gion in the disk where they form, increase with stellar 
mass. Lower mass disks produce failed icy cores, which 
are probably similar to Uranus, Neptune, and the ob- 
served "super-Earths." 

The outward movement of the snow line as stars more 
massive than the Sun reach the main-sequence, and as 
the disk beco mes optically thin, allows the ocean planets 
suggested by iLeger et"aD ([2004) to form in situ. The 
change in disk temperature is only large enough for these 
planets to harbour oceans around stars >2.5M Q . 

Our model includes several poorly determined param- 
eters, which current and future facilities will investigate. 
Whil e there are curren t resolved studies of gaseous disks 
(e.g. lBitner et aLlfeOOTl ). the next generation of telescopes 
such as GMT and ALMA will provide more informa- 
tion on surface density profiles, and how disk properties 
change with stellar mass and age. These studies will 
help to constrain input parameters for our model. We 
have shown that the time dependence of the snow line 
in part determines where gas giant cores form. This re- 
sult should motivate future studies of planet formation 
in disks whose properties change with time. 

The subsequent evolution of isolated cores is beyond 
the scope of this paper, but further work that investigates 
the growth and dynamical evolution of these objects can 
investigate the diversity of resulting system structures. 

Given an initial distribution of disk masses, the prob- 
ability that a star has at least one gas giant increases 
linearly with stellar mass from 0.4 Mq to 3M Q . If the 
frequency of gas giants around solar-mass stars is 6%, 
we predict an occurrence rate of 1% (10%) for 0.4 M Q 
(1.5 Mq) stars. This result is largely insensitive to 
changes in our model parameter s. 

In contrast to the llda fc Linl (|2005| ) model, where it 
is hard to form observable gas giants above 1 Mq , our 
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model predicts a peak at ^3 Mq because we include disk 
and PMS evolution in our snow line derivation. However, 
our model does not include migration, so our prediction 
applies to observable and currently undetectable gas gi- 
ants. 

Though sample numbers are small, it appears that 
observable gas giant frequency increase s with stellar 
mass a cross a wide range of host masses (jJohnson et al.l 
l2007al) . Larger samples of stars that host giant planets, 
particularly low and intermediate-mass stars, will solid- 
ify this result. These studies, and the extension of the 
results to a wider range of semi-major axes, will provide 



a basis for comparison with our model predictions. 
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